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Abstract 

Quadratic variations of Gaussian processes play important role in 
both stochastic analysis and in applications such as estimation of model 
parameters, and for this reason the topic has been extensively studied 
in the literature. In this article we study the problem for general Gaus¬ 
sian processes and we provide sufficient and necessary conditions for 
different types of convergence which include convergence in probability, 
almost sure convergence, L p -convergence as well as convergence in law. 
Furthermore, we study general Gaussian vectors from which different 
interesting cases including first or second order quadratic variations can 
be studied by appropriate choice of the underlying vector. Finally, we 
provide a practical and simple approach to attack the problem which 
simplifies the existing methodology considerably. 

Keywords: Quadratic variations; Gaussian vectors; Gaussian pro¬ 
cesses; Convergence in probability; Strong convergence; Convergence 
in L p ; Central limit theorem 
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1 Introduction 


Quadratic variation of a stochastic process X plays an important role in dif¬ 
ferent applications. For example, the concept is important is one is interested 
to develop stochastic calculus with respect to X. Furthermore, quadratic 
variations can be used to build estimators for the model parameters such as 
self-similarity index or parameter describing long range dependence which 
have important applications in all fields of science such as hydrology, chem¬ 
istry, physics, and finance to simply name a few. For such applications one is 
interested to study the convergence of the quadratic variation. Furthermore, 
a wanted feature is to obtain a central limit theorem which allows one to 
apply statistical tools developed for normal random variables. 


For Gaussian processes the study of quadratic variation goes back to 
Levy who studied standard Brownian motion and showed the almost sure 
convergence 


2 n 

lim Y I ~W±_ 

n —>oo L 2 71 

k =1 



= 1 . 


Later this result was extended to cover more general Gaussian processes in 
Baxter [5] and in Gladyshev |I8j for uniformly divided partitions. General 
subdivisions were studied in Dudley El and Klein & Gine [23] where the 
optimal condition o ( J for the mesh of the partition was obtained in 
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order to obtain almost sure convergence. It is also known that for the stan¬ 
dard Brownian motion the condition o is not only sufficient but also 

necessary. For details on this topic see De La Vega m for construction, and 
m for recent results. Functional central limit theorem for general class of 
Gaussian processes were studied in Perrin m- More recently Kubilius and 
Melichov [25] defined a modified Gladyshev’s estimator and the authors also 
studied the rate of convergence. Norvaisa [29] have extended Gladyshev’s 
theorem to a more general class of Gaussian processes. Finally, we can 
mention a paper by Malukas [25] who extended the results of Norvaisa to 
irregular partitions, and derived sufficient conditions for the mesh in order 
to obtain almost sure convergence. 

The case of fractional Brownian motion with Hurst index H £ (0,1) were 
studied in details by Gyons and Leons m where the authors showed that 
appropriately scaled first order quadratic variation (that is, the one based 
on differences X tk — X tk _ 1 ) converges to a Gaussian limit only if H < |. 
To overcome this problem, a generalisations of quadratic variations were 
used in m, m, and m- The most commonly used generalisation is 
second order quadratic variations based on differences X tk+1 — 2X tfc + X tk _ 1 
which was studied in details in a series of papers by Begyn mm® with 
applications to fractional Brownian sheet and time-space deformed fractional 
Brownian motion. In particular, in [5] the sufficient condition for almost 
sure convergence was studied with non-uniform partitions. The central limit 
theorem and its functional version were studied in [6] and [7] with respect to 
a standard uniform divided partitions. Furthermore, the authors in papers 
mm have studied more general variations assuming that the underlying 
Gaussian process have stationary increments. For another generalisation, 
the localised quadratic variations were introduced in [8] in order to estimate 
the Hurst function of multifractional Gaussian process. These results have 
been generalised in H3E6]. 

The fractional Brownian motion has received a lot of attention in mod¬ 
elling as a (relatively) simple generalisation of a standard Brownian motion. 
However, for some applications the assumption of stationary increments is an 
unwanted feature. For this reason there is a need for extensions of fractional 
Brownian motion. Recent such generalisations are sub-fractional Brownian 
motion depending on one parameter H £ (0,1) introduced by Bojdecki et 
al. [TO], and bifractional Brownian motion depending on two parameters 
H € (0,1) and K £ (0,1] (the case K = 1 corresponding to the fractional 
Brownian motion) introduced by Houdre and Villa [21], and later studied in 
more details by Russo and Tudor [33]. Furthermore, bifractional Brownian 
motion was extended for values H £ (0,1), K £ (1,2) satisfying HK £ (0,1) 
in p]. 
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Recently there has also been interest in general Hermite variations which 
are partially motivated by the contributing paper by Breuer & Major (T2j. 
For results related to fractional Brownian motion, we refer to HH 1301 EH. 
The integrals of fractional Brownian motion were studied in m- Moreover, 
fractional Brownian sheet have received attention at least in papers [421 1431 . 
l39j . Furthermore, generalisations are studied in [38] and [2] . The mentioned 
papers studying general Hermite variations are not focused on practical 
importance of quadratic variations. Furthermore, the proofs are based on 
now well-developed Malliavin calculus. 

This paper aims to take a practical, instructive, and a general approach 
to quadratic variations. Firstly, our aim is to provide easy to check condi¬ 
tions for practitioners which still covers many interesting cases. Secondly, 
with our simplified approach we are able to provide intuitively clear expla¬ 
nations such as the discussion in subsection 13.31 rather than get lost into 
technical details. Finally, to obtain the generality we study sequences of 
general n-dirnensional Gaussian vectors Y n = ^Y^ n \ ..., , where each 

(n) 

component may depend on n, and we study the asymptotic behaviour 
of the vector Y n or its quadratic variation defined as the limit 



k =1 


provided it exists in some sense. As such, different cases such as first or 
second order quadratic variations can be obtained by choosing the vectors 
Y n suitably, and this fact will be illustrated in the present paper. 

We begin by providing sufficient and necessary conditions for the con¬ 
vergence in probability which, applied to some quadratic functional of a 
given process, can be used to construct consistent estimators of the model 
parameters. Furthermore, we show that in this case the convergence holds 
also in L p for any p > 1. We will also apply the well-known Gaussian con¬ 
centration inequality for Hilbert-valued Gaussian random variables which 
provides a simple condition that guarantees the almost sure convergence. 
This condition is applied to quadratic variations of Gaussian processes with 
non-uniform partitions for which we obtain sufficient conditions for the con¬ 
vergence. More importantly, this result is shown to hold for many cases of 
interest and in the particular case of standard Brownian motion, this condi¬ 
tion corresponds to the known sufficient and necessary condition. Compared 
to the existing literature, in many of the mentioned studies the almost sure 
convergence is obtained by the use of Hanson and Wright inequality [2D] to¬ 
gether with some technical computations. In this paper we show how these 
results follow easily from Gaussian concentration phenomena. 
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We will also study central limit theorem in our general setting. We 
begin by providing sufficient and necessary conditions under which appro¬ 
priately scaled quadratic variation converges to a Gaussian limit. To obtain 
this result we apply a powerful fourth moment theorem proved by Nualart 
and Peccati m which, thanks to the recent results by Sottinen and the 
current author [36], can essentially be applied always. We will also show 
how a version of Lindeberg’s central limit theorem for this case follows eas¬ 
ily. Finally, we will use some well-known matrix norm relations to obtain 
surprisingly simple way to obtain a convergence towards a normal random 
variable. More remarkably, it seems that this simple condition is essentially 
the one used in many of the mentioned studies while the result is derived 
case by case. We will also provide a Berry-Esseen type bound that holds 
in our general setting which, to the best of our knowledge, is not present 
in the literature excluding some very special cases. Furthermore, our ap¬ 
proach does not require the knowledge of Malliavin calculus and should be 
applicable to anyone with some background on linear algebra and Gaussian 
vectors. 

To summarise, in this paper we give sufficient and necessary conditions 
for the convergence of guadratic variations of general Gaussian vectors which 
can be used to reproduce and generalise the existing results. Furthermore, 
we give easily checked sufficient conditions how one can obtain the wanted 
convergence results. As such, with our approach we are able to generalise 
the existing results as well as simplify the proofs considerably by relying 
on different techniques. At the best, the methods and results of this paper 
would provide new tools to attack the problem under consideration while 
classically the problem is studied by relying on Hanson and Wright inequality 
together with Lindeberg’s central limit theorem. 

The rest of the paper is organised as follows. In section [2]we study general 
Gaussian vectors and provide our main results. In section [3] we illustrate 
how our results can be used to study quadratic variations. We will consider 
non-uniform sequences and generalise some of the existing results. The main 
emphasis is on first order quadratic variations which is more closely related 
to stochastic calculus while we also illustrate how second order quadratic 
variations can be studied with our approach. We end section [3] with a 
discussion on general quadratic variations which is close in spirit with the 
work by Istas and Lang [22j . Finally, section [4] is devoted to examples. 
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2 Convergence of Gaussian sequences 

2.1 Notation and first results 

Let Y n = (y^ h \ • • • > be a zero mean Gaussian vector, where 

(n) 

each Y£ possibly depends on n. We consider properties of sequences Y n as 
n tends to infinity. Throughout the paper we will also use Landau notation, 
i.e. for a sequences a n and b n we denote; 

• a n = 0(b n ) if sup n > x 1^1 < oo, 

• a n = o (b n ) if lim^oo = 0. 

We also denote a n ~ b n as n —>• oo if lim n _ > . 0O ^ —>• 1. 

On 

We begin with the following definition which is a discrete analogy to the 
similar concepts introduces in ITr . 

Definition 2.1. (i) We say that the sequence Y = (Y n )™ =1 has a 

quadratic variation (Y) if the random variable (Y) := limn^oo Ylk=i (j^k^ 
exists as a limit in probability. 

(ii) The energy of a sequence Y = (Y n )™ =1 is defined as the limit 

n 2 

k=1 


provided the limit exists. 

(iii) We say that Y has 2-planar variation defined as the limit 

2 


n m 

T(y) := lim VVfe 

n,m— Kx> [_ 

k=l j =1 


yWyW 
k 3 


provided the limit exists. 
We will also denote 




k =1 


Y^Y-ElYf ’' 2 


( 2 . 1 ) 


for the centered quadratic variation. 
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Example 2.1. Let X = {X t ) te \n T \ be a centred Gaussian process, ^ 
and let A^,X denote some difference of the Gaussian process X. By setting 


(n) _ A k X 


_ 

^ h - 


nLxy we have E ( Y k 


( n )\ _ ^,-1 


= n 


and 


K = 


A k X 


1 n 

n^E(A fc X) ; 


- 1. 


Hence by setting A^.X = X*n — X^ we obtain the first order quadratic 
variation with respect to uniform partition. Similarly, by setting A&X = 
Xfi ^ — 2Xt" + X^ i we obtain the second order quadratic variation with 
respect to uniform partition. 

The Euclidean distance of the vector Y n = (Yj n \Y^ n \ ..., Yn 1 ' 1 ) is given 
by 


N 




k =1 




Note that the norm || • ||2 also depends on the dimension n which we will 
omit on the notation. We will denote by r( n ' the covariance matrix of the 
vector Y n , i.e. the n x n-matrix with elements 


i(”) 


T = E ( Yr'Y,. 


'(n) v (n) 


3 


( 2 . 2 ) 


Note now that with this notation the energy of a process Y is simply the limit 
of the trace of the matrix T, i.e. e(Y) = lim^-^ trace(T n ). Similarly, Y n 
has a quadratic variation if the limit ||y n, ||| converges as n tends to infinity. 
Recall next that the Frobenius norm of a matrix is 

given by 


||r (n) || F = 


v 


E 

i,j =1 


n(n) 

ij 


E 


yMy'l 71 ) 

k j 


Hence we have 

n n 

lim ||r (n) ||| = lim \. 

n—> 00 n—too 2 —' L 

k =1 j =1 

We will later show that in interesting cases we also have 

n m ^ n n 


11 2 


(2.3) 


lim 

n,m—> 00 


EE[ e 


yWyH 
k J 


= lim 

n—>-00 


EE[ e 


y( n )y( n ) 

^ j 


11 2 


/c=i j=i fc=i 1=1 

which, in view of (12.31) . shows that the 2-planar variation T(E) is given by 


T(y) = lim ||r( n )| 


F- 
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The following first result concerns convergence in probability. The proof 
follows essentially the arguments presented in [05] and is based on cumulant 
formulas for Gaussian random variables. The main difference is that we 
also prove the convergence in L p for any p > 1 which have some important 
consequences in stochastic analysis (to be detailed in a forthcoming paper) 
while in [05] the authors considered only L 2 -convergence. The proof follows 
the ideas presented in [05] but we will present the key points for the sake of 
completeness. 

Theorem 2.1. Let Y = (Y n )^ =1 be a sequence of Gaussian vectors with 
finite energy. Then quadratic variation exists as a limit in probability for 
every t > 0 if and only if the sequence (Y n )^ =1 has 2-planar variation. In 
this case, the convergence holds also in L p for any number p > 1. 


Proof. Denote 


= iiml = E T« 


dn) 


k =1 


We start by showing that 


Z p = 


E(n 


k= 1 


d n ) 


is uniformly integrable for any p > 1. By Minkowski’s inequality for mea¬ 
sures we get 


d”) 


E E ( y > 


Kk=1 


< 


n 

E E (h” ) ) 


2 P 


k =1 L 


c„E E n 


r(n) 


by the fact that Y^ is Gaussian. Now this upper bound converges to e(Y), 
and consequently the quantity ( Y^-=\ ) is uniformly integrable for 


any p > 1. Now we have 

E(Z n ~~ Z m y = EZ 2 + E Zf n — 2E (Z n Z m ), 

where 

n m 

^{ZnZm) = E y± k 

k= 1 j=l 
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Here we have 


E 




= 2 


E 


( y Wy( m ) 
\ k j 


2 + E^ n) ) 2 E(h} (m) ) 2 . (2.4) 


Hence we have 


n m 


EE e 



Consequently, 


= 2 


n m 

EEK 


y( n )y( m ) 


fc=l j=l 


1 2 


+E E 

fc=i 


y 


(n) 


: E E 

k =1 


y 


(n) 


(2.5) 




2 EE[ E (^ , h w ) 


fc=l j=l 


n n 

+E E (y fc (n) ) 2 x^E(y fc (n) ) 2 , (2.6) 

k =1 /c=l 


where 

n n 

X>(y fc (n) ) 2 x ^E(y , (n) ) 2 -»• e(yf 

k =1 

by the fact that Y has finite energy. Recall also that 2-planar variation is 
given by 

n m 2 

T < y >=„“!5oo E E [ E ■ 

k= 1 i=i 

Hence relation (12.61) implies the result. Indeed, assuming that Z n converges 
in probability, then uniform integrability implies that (Y) € L p and the 
convergence holds also in L p . Hence E (Z n Z m ) converges, and (12.61) implies 
that 2-planar variation exists. Conversely, if 2-planar variation exists, then 
(12.61) implies that E (Z n Z m ) converges to the same limit as E (Z%) which 
concludes the proof. □ 


Remark 2.1. It is straightforward to check that the IJ 1 -convergence takes 
place also in a continuous setting of Russo and Vallois |45j . 


The following theorem gives condition when the quadratic variation is 
deterministic. It seems that this is indeed true in many cases of interest. 

Proposition 2.1. Let Y = (Y n )^ =1 be a sequence of centred Gaussian 
vectors such that Y has finite energy. Then quadratic variation exists as a 
limit in probability and is deterministic if and only if the 2-planar variation is 
zero. In this case (Y) equals to the energy of the process and the convergence 
holds also in L p for any p > 1. 


Proof. The result follows directly from 

lim E (Z n - e(y )) 2 = E (Z%) - e(Y) 2 
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and the relation (12.61) with m = n. Finally, the convergence in L p follows 
directly from Theorem 12.11 □ 


Remark 2.2. A generalisation of quadratic variation is a-variation which 


is defined as a limit of Ylk=i 


Y, 


(n) “ 


say that Y has finite a-energy if the limit 


provided the limit exists. Similarly, we 


*cn := lim 

n .—Vnn z J 


(n) 


k =1 


exists. It is straightforward to show that if the sequence Y 
finite a-energy and 



k =l 


(W=i has 


converges in probability, then the convergence holds also in L p for any p > 1. 
In this case however, the concept of 2-planar variation becomes much more 
complicated. 


2.2 Almost sure convergence 

In this subsection we address the question when the convergence takes place 
almost surely. The key ingredient for our results is the concentration in¬ 
equality for Gaussian measures and by use of this inequality we show that, 
rather surprisingly, the quadratic variation always converges to the energy 
of the process whether or not the energy is finite provided that 2-planar 
variation vanishes. 

Before stating our results we recall the following concentration inequality 
taken from [JJ for Gaussian processes. We present the result using our 
notation. 

Lemma 2.1. Let Y n be a Gaussian vector with covariance matrix , and 
denote T = yjtrace (T( n )). Then for any h > 0 we have 

p(ini 2 -ri>/ l )< rap (- ij Aj 1 -). (2.7) 

Remark 2.3. The result holds for any Hilbert-valued Gaussian random 
variables, not only finite dimensional. Similarly, for general a-variations one 
can use the concentration inequality 

P(|||X|| B -E||A:||B|>/i)<exp 


10 








where (. B , || • ||g) is a Banach space, X is a Banach-valued Gaussian random 
variable and a 2 = sup^gg/.MMi^ E,L(X) 2 . Applying this to W 1 equipped 
with the norm || • || Q together with Riesz representation theorem one has 

p(|||y|| a -E||yy > h) < ex P , (2-8) 

where u 2 = sup^H^ ELi £"= 1 afcajE[y fe ^] with I + | = 1. 

We now turn to address the question when one obtains almost sure 
convergence. Clearly, the idea is to find an upper bound on || r ( n )|| 2 , say, 
|| r ( n )|| 2 < (j>(n). Then one can hope that the bound is good enough such 
that 

°° / e 2 \ 

(29) 

from which we obtain the almost sure convergence immediately by Borel- 
Cantelli Lemma. 


Theorem 2.2. LetY 
||r( n )|| 2 -> 0 and 


(Y n )° < L 1 be a sequence of Gaussian vector such that 



( 2 . 10 ) 


Then, as n —>• 00 , 



E E ( r i= (n) ) 


2 


-)• 0 


( 2 . 11 ) 


in probability and in L p for any p > 1. Furthermore, the convergence holds 
almost surely for any partition satisfying ||r <n )||2 = o (j-i— ). 


Proof. The convergence From Lemma 12.II we get that 


\ 

M* 

111 

to 

1 

n 9 

L E (d n) V 

\ 

II 

?r* 

II 


in probability follows immediately from Lemma [2.11 and the almost sure con¬ 
vergence follows by applying Borel-Cantelli Lemma. Now the convergence 
112.111) follows from decomposition \a—b\ = \y/a — \/b\(y/a+Vb) together with 
the uniform integrability condition (12.101) which also implies convergence in 
L p . □ 
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Remark 2.4. Note that if ||r( n )|| 2 —> 0, we have 


\ 

M* 

IE 

to 

1 

-•> 

n 9 

\ 

k =1 \ 

II 


in probability even if sup n>1 Ylk=i ® {^k"^ = °°> be. we h ave convergence 
in probability (or almost surely) while the energy might be infinite. For 
example, it can be shown that this is the case for fractional Brownian motion 
B h with Hurst index H G (0, and its non-scaled quadratic variation, i.e. 
the one corresponding to a vector = B[^ — B^_ . 

Remark 2.5. Note that for finite energy processes such that 2-planar vari¬ 
ation vanishes this result shows that 


lim 

n,m—Kx) 


n rn 

k= 1 j =1 


E 


yWyW 
k j 


lim VV E 

n—>oo l—s Z—' L 

k =1 j =1 


k j 


2 


or more compactly, Y(Y) = lim^-^x, ||r^ n ^||jr. Indeed, from the well- 
known relation ||A|| 2 < ||^4 ||f we obtain that if lim n _ > , 00 ||rl n l—>• 0, 
then ||r ( n )|| 2 —> 0. Consequently, Theorem 12.11 implies that 2-planar vari¬ 
ation vanishes. This answers to question raised in [451 Remark 3.12] in 
our discrete setting. Similarly, one can use general concentration inequality 
(12.81) to give analogous answer in a continuous case. 


To compute the spectral norm || 2 , or equivalently the largest eigen¬ 
value, can be a challenging task. One way to overcome this challenge is to 
use Frobenius norm || • \\p which provides an upper bound and is easier to 
analyse. Unfortunately however, it provides quite rough estimates even in 
a simple case of standard Brownian motion as will be shown in subsection 
an a way to obtain general conditions is to use matrix norm || • ||i which 
is also the main approach applied in the literature. This is the topic of the 
next general theorem. The proof is based on some well-known relations for 
matrix norm and we do not claim originality here. 


Theorem 2.3. Let Y = (Y n )ff =1 be a sequence of Gaussian vectors such 
that (12.101) holds. Furthermore, assume there exists a function cj)(n) such 
that 


max 

3 


n 

y( n )y( n ) 
k j 

k =1 


E E 


< (f{n). 


If 4>(n) —>• 0, then the convergence 



n 


-> 0 
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holds in probability and in L p for any p > 1. Furthermore, if 

m = ° (ii) ■ 

Then the convergence holds almost surely. 


Proof. Recall the well-known bound for matrix norm ||H ||2 given by ||H ||2 < 
y/WAUAlU where ||^4||i = maxi< 3 -< n Xfc=i Wj\ and \\A\\^ = maxi<j< n Xfc=i \ a jk\- 
Hence symmetry of gives ||rW || 2 < ||r( n )||i, where 


||rWib = 


i = max 


Ei E 


y( n )y( n ) 
k j 


k =1 


This proves the claim together with Theorem 12.21 


□ 


The following final result of this section can be useful for stationary 
sequences. 

Theorem 2.4. Let Y = (Y n )™ =1 be a sequence of Gaussian vectors such that 
(12.101) holds. Moreover, assume there exists a positive symmetric function r 
such that 


E 


y( n )y( n ) 
k j 


<r(k-j), k,j = 1, 


,n 


and assume that there exists a function (j>(n) such that 

n 

r(k) < 4>(n). 


k =i 


If <f>{n) —>• 0 as n tends to infinity, then the convergence 


k =1 k =1 


(n) 


1 2 


holds in probability and in L p . Furthermore, the convergence holds almost 
surely provided that 4>(n) = o ( J . 


Proof. Now for any j > 1 we have 


EI e 


y( n )y( n ) 
k j 


3-1 


n-j 


<J2 r (k-j) + r{k) < 2 E--(t) 


k =1 k =1 

from which the result follows. 


k =o 


fc=i 


k =0 


□ 
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2.3 Central limit theorem 


In this section we provide sufficient and necessary condition for the central 
limit theorem (CLT) to hold. More importantly, we obtain two simple corol¬ 
laries where the first one gives a version of Lindeberg’s CLT for quadratic 
variations, and the second one gives simple to check condition which actu¬ 
ally holds in most of the studies cited in the introduction. In particular, 
the second corollary can be used to simplify the used techniques (namely, 
the Lindeberg’s CLT) considerably. Our necessary and sufficient condition 
is based on the fourth moment theorem by Nualart and Peccati m- Hence 
we begin by recalling some basic facts on Wiener chaos. For more details 
we refer to monographs [351 Sol EH- 

Let now X be a separable Gaussian centered process and let q > 1 be 
fixed. The symbol TL q denotes the qth Wiener chaos of X, defined as the 
closed linear subspace of L 2 (fl) generated by the family {H q (X(h)) : h € 
n,\\h\\ H = 1}, where H q is the qth Hermite polynomial. The mapping 
I*{h® q ) = H q (X(h)) can be extended to a linear isometry between the 
symmetric tensor product TL Qq and the gth Wiener chaos TL q , and for any 
h 6 T-L Qq the random variable I^{h) is called a multiple Wiener integral of 
order q. 

Remark 2.6. If X = IF is a standard Brownian motion, then 7~L is simply 
the space L 2 ([0, T], df). In this case the random variable I' Y (h) coincides 
with the g-fold multiple Wiener-Ito integral of h (see [35] 1. 

Remark 2.7. Let X be a separable Gaussian process on [0, T\. It was 
proved in [46] that X admits a Fredholm integral representation 

X t = [ K{t,s)dW s , 

Jo 

where K € L 2 ([0,T] 2 ) and IF is a Brownian motion, if and only if 
J 0 T EX 2 du < oo. Furthermore, it was shown that this representation can 
be extended to a transfer principle which can be used to develop stochastic 
calculus with respect to X. In particular, the transfer principle can be 
used to define multiple Wiener integrals as multiple Wiener integrals with 
respect to a standard Brownian motion. This definition coincides with the 
one defined via Hermite polynomials. 

Finally, we are ready to recall the following characterisation of conver¬ 
gence towards a Gaussian limit. 

Theorem 2.5. J5?|/ Let {F n } n >i be a sequence of random variables in the 
qth Wiener chaos, q > 2, such that lim^^ooE(F 2 ) = a 2 . Then, as n —» oo, 
the following asymptotic statements are equivalent: 
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(i) F n converges in law to W(0,a 2 ). 

(ii) E Ft converges to 3a 4 . 

Remark 2.8. In this paper we are studying quadratic variations of Gaussian 
sequences. Hence, thanks to Fredholm representation from [46], such objects 
can be viewed as a sequences in the second chaos. 

Remark 2.9. The case of the second chaos was studied in details in [32,23] 
where the authors characterised all possible limiting laws. More precisely, 
authors in [32] proved that if a sequence in the second chaos converges in 
law to some random variable F, then F can be viewed as a sum of normal 
random variable and an independent random variable living in the second 
chaos. 

Remark 2.10. It was proved in [36] that instead of fourth moment one can 
also study the convergence of ||.D.F n ||^ in L 2 , where D stands for Malliavin 
derivative operator, and this approach have turned out to be very useful in 
some cases. For our purposes however it seems that working with the fourth 
moment is more convenient. 


Finally, we recall the following Berry-Esseen type estimate taken from 
[301 Theorem 11.4.3]. 

Theorem 2.6. Let {F n } n >i be a sequence of elements in the qth Wiener 
chaos such that E( F,' 2 ) = 1 and let Z denote a standard normal random 
variable. Then there exists a constant C q depending only on q such that 


sup 

irER 


P(F n < x) — P(Z < x) < C q ^EF^ — 3. 


We are now ready to prove our results. We begin with the following 
auxiliary technical lemma whose proof is postponed to the appendix. 

Lemma 2.2. For V n given by m we have 


EV 2 = 2 £ (E 

k,j =1 


^( n )^( n ) 
k j 


and 


EV* = 12 


-I 2 


E( e 

k,j =1 


yd n )yd Tl ) 

k j 


+ 24 E 

i,j,k,l =1 


^( n )^( n ) 
k j 

E 

IT 

Jr 

IT 

E 

5 

J 

E 

i i 

-"I 

i_i 
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With the help of this lemma we are ready to proof our main theorem. 

Theorem 2.7. LetY = (Y n )ff =1 be a sequence of Gaussian vectors such that 
for every n > 1 the elements Yf , k = 1,..., n belong to the first Wiener 
chaos, and let denote the covariance matrix of Y n with eigenvalues 
A". Xlf ,..., A”. Then, as n tends to infinity, the following are equivalent. 


Vn 


y/Var(V n ) 


—> m, i)? 


E E 

i,j,k,l =1 


y( n ) Y'fa) 

A; j 

E 

5 

2 

E 

-"I 

E 

& 
i_i 


= o ( Var{V n ) 2 ) , 


E< A 

k =1 


4 _ 

k) ~ 0 


1 2\ 


Ymf 

Lk=l 


Proof. By assumption we are able to use fourth moment theorem 12.51 from 
which the equivalence of items (1) and (2) follows with the help of Lemma 

12.21 To obtain equivalence of (1) and (3), it is well-known that Ylk =l 
can be decomposed as 


Y, 


(n) 


n rx n 

EM 


k =1 


i i 


k =1 


(n) 

k 


where A^ are the eigenvalues of the covariance matrix 1^”) and are 
independent standard normal random variables. Denoting 


\ r 
* i 


(n) 


u = E 

k =1 

and using Lemma 12.21 again we obtain 

n 

ED n 2 = 2^ 

and 


_ V 

Ai 


ni2 


k =1 


ED n 4 = 12 

which concludes the proof. 


Ei A 

k =1 


ni2 

fcJ 


+ 6^[A 


n ]4 

fcj 


fc=l 


□ 
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As a simple corollary we obtain the following result which corresponds 
to Lindeberg’s CLT and is mainly used in the references given in the intro¬ 
duction. 


Corollary 2.1. Assume that assumptions of Theorem \ 2. 7| prevail. If 


A*(n) := max |A£| = o 
k=l,...,n 


n —f oo, 


then 


V n 


VVar(V n ) 


Af(0,l). 


Proof. We have 

n n 

k =1 ’"' ,n k =1 

and since Car(IL n ) = f>2k=i(^k) 2 > the result follows at once. □ 

Remark 2.11. Note that since Lindeberg’s CLT can be proved without the 
theory of Wiener chaos, the above result is valid for arbitrary sequences Y n . 


Finally, the following theorem justifies that in many cases it is sufficient 
to find an upper bound for A*(n), or even for maxi<j< n E Yj n> . 

While the proof follows from simple relations for matrix norm, the result turn 
out to be very useful for many practical applications. In particular, the fol¬ 
lowing result covers many of the cases studied in the literature. Furthermore, 
in this case it easy to give a Berry-Esseen bound. 


Theorem 2.8. Let the assumptions of Theorem 1 2. 7| prevail, and assume 
that Y n is a Gaussian vector with finite non-zero energy. Then there exists 
a constant C > 0 such that 


sup 


p 


v n 


< x — P (Z < x) 


VVar(V n ) 

where Z is a standard normal random variable. Hence if 


< Cy/n\*{n). 


then 


v n 


y/Var(Vn) 


then 


Vn 


y/Var(V n ) 


A-(„)=„(_!, 


Af(0, 1). In particular, if 

n 

y |e 


max 

1 <j <n 


y^( n )y^( n ) 
k j 


k =1 


= O 


Af(0,l). 
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Proof. Recall the trace norm is given by ||r( n l||* = Y^k=i^k- By Cauchy- 
Schwartz inequality we get the well-known matrix norm inequality 


n 



By assumption, Y n has finite non-zero energy. Hence by observing that 
lim n ^ 00 ||r( n )||* = e(Y) > 0, we obtain that for large enough n we have 


o < c< V^||r( n )|| F . 


Next we observe that y/Var(V n ) = ||Tl”)||i?. Now the fourth moment of 

1 s is given by 
y/Var(V n ) 

EV* 3 | 3ELiM] 4 

(EK ?) 2 2 (£LiMl 2 ) 2 


so that 


Hence the Berry-Esseen bound follows from Theorem 12.6 


EK 4 3 < 3[A*(n )] 2 < 3 2 

(EI 4 2 ) 2 2ER n 2 2 


□ 


Remark 2.12. Note that the convergence towards normal random variable 
follow also from Corollary 12.11 which does not rely on the theory of Wiener 
chaos. However, for a sequence living in the second chaos we also obtain a 
Berry-Esseen bound. 


3 Application to quadratic variations 

In this section we apply the results to quadratic variations of Gaussian 
processes. We begin by giving a general results for generalised variations 
which will be illustrated in the particular cases of first and second order 
quadratic variations. 

We consider arbitrary sequences of partitions vr n = {0 = fg < tf < 
... < = T}, where N(7r n ) denotes the number of points in the 

partition. For notational simplicity, we drop the super-index n and simply 
use tk instead of The mesh of the partition is denoted by | 7 r n | = maxjffc — 
tk-i,t k € 7 r n /{ 0 }}. We also use m( 7 r n ) = min {4 - 4-1,4 € 7 r n /{ 0 }}. 
Throughout this section we assume that 

^ T < fc(knl), n> 1 (3.1) 

m(7T n ) 

for some function k. Obviously, usually the partition is chosen such that 
fc(|7T n |) < C < oo. 
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3.1 First order quadratic variations 


In this subsection we apply the results of previous section to study first order 
quadratic variations of Gaussian processes which is our main interest due 
to its connection to stochastic analysis. Throughout this subsection we also 
use the metric defined by the incremental variance of X, i.e. 

d x (t,s) = E(X t -X s ) 2 . 


Definition 3.1. Let X = (X t ) te ^ T ] be a centred Gaussian process. We say 
that X has first order ^-quadratic variation along n n if 

J7 , -^tk- 1 )^ 

Vi(n n , (f>):=2_^ 


t?TTn 


tk— l) 


(3.2) 


converges in probability as |7r n | tends to zero. 

Remark 3.1. A natural choice for the function ([> is such that 


lim 

\^ n \—^0 


AT(7r n )-l 

E 

k=1 


E (X tk -X tk _^ 
J'k ^k—l) 


= K < oo. 


(3.3) 


In particular, in many interesting cases one has dx(t , s ) ~ r(t—s) as \t— s| —> 
0 for some function r. In this case a natural choice is 4>(x) = 

To simplify the notation we denote V\ (n n , <f) = V\{tt n ,4>) — EVi(7r n ,<f>). 
We also use A tj = tj — tj- 1 . We will begin by giving the following general 
theorem which generalises main results of [28] by allowing us to drop some 
technical assumptions. The result follows directly by uniting and rewriting 
Theorems 12.21 and 12.81 for sequence Y 7 n] = X .‘ k 

VWfc-*fc-1) 

Theorem 3.1. Let X be a Gaussian process. Assume that 

N( 7Tn)-l 


max 


E 


1 


i<3<jv(7r„)-i v /0(A4)^(Ai :? -) 

for some function iL(|7r n |). 


|E[(W jt -W fc _ 1 )(W,-W j _ 1 )]| < H (|7T„|) 

(3.4) 


(i) If H(\ir n \) —> 0 as |7r n | tends to zero, then convergence 
N M-i (x _ x )2 N ^~ 1 E ( X t , - X tk .) 2 

\ A \ L k L k — 1/ \ A \ L k L k — 1/ v ft (q 

^ m-tk-i) ^ m-tk-i) 

holds in probability. If H( |vr n |) = then convergence P 

holds almost surely. In these cases the convergence holds also in L p 
for any p > 1 provided that P holds. 
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(ii) Furthermore, assume that 


lim M V lE <^V =g > o. 

Kn|—>0 *P\tk tft — 1) 


Then there exists a constant C > 0 such that 

Vl (7T n ) 


(3.6) 


sup 

ieER 


p 


Var(Vi(n n )) 


< x — P(Z < x ) 


<Cy/N(Tr n )H(\n n \), 


where Z is a standard normal random variable. In particular, if 
H(\ir n \) = o (^N(ir n )~^, then 


Vl (7T n ) 


V ar(Vi(ir n )) 


AA(0,1). 


Remark 3.2. In [28] the author studied a particular class of Gaussian pro¬ 
cesses while here we consider arbitrary Gaussian process. Similarly, in [28] 
the main result was derived by using some technical computations under as¬ 
sumption (13.411 together with several additional technical assumptions. Here 
we have shown that (13.41) is the only needed condition which generalises the 
class of processes considerably. Similarly, we have been able to simplify the 
proof since we have shown that such results follows essentially from Gaussian 
concentration together with some matrix algebra. 


Next we will provide some sufficient conditions which are easy to check. 
A particularly interesting case for us is Gaussian processes such that the 
function 

d(s,t) =P(X t -X s ) 2 

is C 1,1 outside diagonal. Note that a sufficient condition for this is that 
the variance E2f t 2 is C l and the covariance R of X is C 1,1 outside diagonal. 
Furthermore, note that this assumption is satisfied for many cases of interest. 

The first theorem gives a general result in a case of bounded derivative. 

Theorem 3.2. Let X be a continuous Gaussian such that the function 
d(s,t) = E(Xt — X s ) 2 is C 1,1 outside diagonal. Furthermore, assume that 
there exists a positive function f(s,t) such that 

|dstd(M)l < f(s,t) 


and 


pv 

SU P / f(s,t)dt<oo. 
s,ve[o,T] Jo 
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Assume there exists a function H(\ir n \) such that 


dfe.ij-l) + Atj 

- £ HiM) - 

Then the result of Theorem \S. 1\ holds with function H(\n n \). 

Proof. Now for j k we have 


E[(X tk -X tk _ 1 )(X tj -X tj _ 1 )\ = 


ftk rtj 
i)}= / d st d(s,t)dsdt. 

Jtk-i Jtj —i 


Hence we have 


E F 


^( n )yd n ) 

k J 


k =1 


d(tj,tj— i) 


1 iV(7r n )-l 


0 (Afj) y/^Atj) fe= ^- y/<j>(At k ) 


E 


rtj 

/ / |<9 st d(s,f)| dsdt 

Jtk-i Jtj-1 


< d{tj,tj_ i) 


JV(ir„)-l 


E 




f(Atj) cfiAtj) vWA T k 


= f L f 3 \d st d(s,t)\dsdt 

cj Jtk-1 Jtj -1 


< 


-i/ <p(Atj) 

l g7*‘ f Mm) | 

^ At ;> .=«#/ V. v. 


<MAfj) 
Now here 


dsdf. 


JV ( 7r ™)- 1 Pt fc Pt,- 


^2 [ f \d st d(s,t)\dsdt 

= 1 ,fc^j J^-l Jtj-1 

f ( \d st d(s,t)\dsdt+ j f \d st d(s,t)\dsdt 


= O (Afj) 

by Tonelli’s theorem and assumptions. The claim follows at once. □ 

Remark 3.3. The convergence depends now on sup 1<fcj<J Y( 7rri )_i j • 

However, in any natural chosen sequence partition we have sup n>1 A;( 7 r n ) < 
oo which guarantees sup 1<fej - <JV ( 7rn )_ 1 j < oo for many functions cj>. 

For example, this is obviously true for power functions <£>(a;) = x 7 which is 
a natural choice in many cases. 
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As an immediate corollary we obtain the following which again seems 
to generalise the existing results in the literature. Most importantly, the 
following result shows how the lower bound for the variance plays a funda¬ 
mental role. Furthermore, a standard assumption in the literature is that 
d(t,s) ~ \t — s | 7 for some number 7 € ( 0 , 2 ) which in particular covers the 
fractional Brownian motion and related processes. In the following the struc¬ 
ture of the variance can be more complicated. For simplicity we will only 
present the result in the case of bounded function k(ir n ) while the general 
case follows similarly. 


Corollary 3.1. Let the notation and assumptions above prevail. Further¬ 
more, assume that there exists a function r such that d(t , s) ~ r(t — s) as 

!*-«!->• 0. Let sup re >! SUp 1 < feJ < JV ( 7rn )_ 1 ^1^) < 00 > SU Pn>l fc(TTn) < OO, 

and put cj)(Atj) = '^ . 

(i) If (\n n \) 2 = o(r(\ir n \)), then 

\Vi(TT n , (f)\ ->• 0 


in probability and in L p for any p > 1. The convergence holds almost 
surely for any sequence = o ( 1 ^ 7 ) • 

(ii) There is a constant C > 0 such that 


sup 

ccEM 


Vi (n n ) 


V ar(\f(7r n ')) 



P (Z < x ) 


< C 



Kknl)’ 


where Z is a standard normal random variable. In particular, if 

3 

(KnI) 5 = o(r(|7r n |)) ; then 


, ^AA(0,1). 

V ar(Vi(7r n )) 


Proof. The claim follows immediately from Theorem 13.21 bv noting that our 
choice of function <f> guarantees condition (13.611 . □ 


We end this section with a following result that recovers the case of 
fractional Brownian motion and related processes. Note again that our 
technical assumptions are quite minimal. 

Theorem 3.3. Let X be a continuous Gaussian such that the function 
d(s,t) = E (Xt — X s ) 2 is in C 1,1 outside diagonal. Furthermore, assume 
that 

\d st d(s,t)\ = O (\t- s\ 2H ~ 2 ) 
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for some H £ (0,1), H / \ and assume there exists a function H( | 7 r n |) such 
that 


max 

l<j<N{lTri)-l 


d(tj,tj- 1 ) + [Atj] min ( 1,2f/ ) sup 1 < feJ < JV(7r7i) _ 1 
4>(A tj) 


V<^(Atfc) 


< H (17T n |). 


Then the result of Theorem, \S. 1\ holds with function H(\ir n \). 


Proof. Note that the case H > A follows directly from Theorem 13.21 Let 
now H < b. Using Fubini’s Theorem we have 


|O s td(s, t )| dsdt 


< C 


T ( 

J tj —l 


(T - s) 2H ~ 1 ds + C I (tj - s ) 


2H-1 


ds 


'tj -1 




for some unimportant constants C which vary from line to line. Here we 
have used the fact that for positive numbers a, b and 7 € ( 0 , 1 ) we have 
|a 7 — 6 7 | < \a — by. Treating integral J ^ -1 1 \d s td(s, t)\ dsdt similarly the 

result follows. □ 


Remark 3.4. We remark that the case H = \ can be treated similarly. In 
this case one obtains a condition 


max 

l<j<N{TT n )-l 


d(tj,tj- 1 ) + Atj\ log Atj\ sup^^jv^).! 




<t>(A tj) 


< H(\n n |). 


Remark 3.5. It is straightforward to give a version of Corollary 13.II in this 
case also. Indeed, the only difference is slightly different exponents in the 
case H < 


3.2 Second order quadratic variations 

In this subsection we briefly study second order quadratic variations. In 
particular, we reproduce and generalise the results presented in papers [5] 
and [ 6 j. We will present our results in slightly different form. However, 
comparison is provided in remark f3.61 

Usually second order quadratic variation on [0,1] is defined as the limit 
of £ 2 =i (x fc+i — 2 Xk + X k-i ). To generalise for irregular subdivisions 

\ n n n J 

Begyn introduced and motivated [5j a second order differences along a se¬ 
quence 7r n as 

AX k = Atk+iX ti ,_ 1 + A tkX tk+1 — (A t k+ i + A t k ) X tk , 
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and we study the second order quadratic variation defined as the limit 


V2(n n ) 


N y ’ 1 At fc+ i (AI fc ) 2 


Again we use short notation 

V2(7r n ) = V 2 (n n )-EV 2 (7T n ). 

We also assume that the derivative g f], d 2 t R{s, t) of the covariance function 
R of X exists outside diagonal and satisfies 


d 4 

d 2 sd 2 t 


R(s , t) 


C 

~ \t — s| 2+ ')' 


(3.7) 


for some number 7 € (0,2). Finally, we make the simplifying assumption 
sup„ k(7r n ) < 00 on the function k. Hence it is also natural to assume 


E(AX fc ) 2 

sup-^ 

j,k E (AX k ) 2 


< 00 . 


(3.8) 


In particular, the assumptions made in [5] implied 

E(Aljt ) 2 ~ (A t k+ i) 2 (A t k ) 2 (At k+1 + At k ) (3.9) 


in which case (13.81) is clearly satisfied. 

Theorem 3.4. Let all the notation and assumptions above prevail. 
(i) If H{ 7T n ) := max 1 <j<jv( 7r „)-i conver 9 es to zero > then 


0 


V 2 (n n ) 

in probability and in L p for any p > 1. Furthermore, the convergence 
holds almost surely provided that H( 717 ) = o 

(ii) There is a constant C > 0 such that 


sup 

fcEK 


v 2 K 


Var(V 2 (7r n )) 


< x I — P (Z < x ) 


< Cy/N(ir n )H(\n n \), 


where Z is a standard normal random variable. In particular, if 
H{TT n ) = o f[ 7 V( 7 r n )] -2 Y then 


V 2 {n n 


Var V 2 (7r, 
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Proof. Denote 


Y(n) _ 

yjE(AX k ) 2 

For some constant C we have 


JV(ir„)-l 

E 

k =1 


E 


y( Tl )y( n ) 
/c J 


I | AT(7r n )-l 

£ |E[A *‘ Ajyi 


by (13.81) and boundedness of k(Tr n ). Furthermore, it was proved in [5] that 
boundedness of fc(7r n ) together with (13.71) yields 


N(n„)-1 

J2 \^[^X k AX j }\<C\n n \ 4 -' y 

k=1 

for a constant C. This gives a bound for 

iV(7r„)-l 

E | E 

k=\ 

from which the result follows immediately by combining Theorems 12.21 and 
f278l □ 

Remark 3.6. To compare our result with the ones provided in papers |5] 
and [6], first note that we were able to reproduce and generalise the main 
theorem of [5] although we gave our result in a slightly different form. Indeed, 
in [5J several additional technical conditions were assumed to ensure the 
asymptotic relation (13.91) while here we have worked with general variance. 
This is helpful since the message of our result is that basically one has to 
only check the asymptotic behaviour of the variance, and (13.71) quarantees 
the upper bound for the corresponding matrix norm. Furthermore, the 
central limit theorem in [6] was proved only for uniformly divided partition, 
and the central limit theorem was proved under more restrictive technical 
conditions, similar to those in [5], and by finding a lower bound for the 
variance in order to apply Lindeberg’s CLT. Here we have proved that such 
result holds also non-uniform partitions and the result follows easily from 
the computations presented in |S] together with Theorem 12.81 Finally, here 
we also obtained Berry-Esseen bound. In particular, under (13.91) we obtain 
bound proportional to y|7r£|. 


y ( n )y( n ) 
k j 


3.3 Remarks on generalised quadratic variations 

We end this section by giving some remarks on generalised quadratic varia¬ 
tions introduced by Istas and Lang [22]. 
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Let now a = (ao,a±,... ,a p ) be a vector such that J2k=o ak = where 
p is a fixed integer. Let also A be a fixed small number and consider time 
points tf. = kA, k = 1,..., n. The a-differences of X is given by 

p 

AaX J = a ' kXt 'j+k > j = 0,l,...,n-p. 
k =0 

In [22] the authors considered stationary or stationary increment Gaussian 
processes and studied generalised a-variations defined as a limit of 

, 1 ^ (A a X k ) 2 

Now since X is either stationary or has stationary increments, the function 
d(t,s) depends only on the difference \t — s|. The assumption in [22] was 
that the function v(t) = d(0,t) is 2D times differentiable ( D is the greatest 
such integer, possibly 0), and for some number 7 € (0, 2 ) and some constant 
C > 0 we have 

v ( - 2D \t) = v ( - 2D \0)+CD + r(t), 

where the remainder r satisfies r(t) = o(t 7 ). The main results in [22] was 
that under different set of assumptions one can obtain Gaussian limit with 
some rate with a suitable choice of the vector a, although in some cases one 
has to assume the observation window nA to increase to infinity. Obviously, 
by using the result of this paper we could reproduce and generalise, at 
least to cover more general variances as in here for first and second kind 
quadratic variations, these results together with a much simplified proofs. 
Instead of getting lost into technical details we wish to give some remarks 
and explanations. In [22] the main message was roughly that larger the value 
of D and s, then larger one has to choose the value p, i.e. taking account 
more refined discretisation. However, the reason for this can, and at least 
was for the present author, be lost in the technicalities. In a nutshell, the 
idea is to try to find a discretisation vector a so that 


1 


max-- - = = 

1 <j<n n 2 ^~E[A a Xj 


2 


n—p 

£ 


k =1 


1 

^E[A a X k } 2 


\E[(A a X k )(A a X j )]\=ofn-^ 


from which we obtain almost sure convergence and central limit theorem 
(with \JVar (V (a, n) — EH (a, n)) as normalisation so one is left to analyse 
the asymptotic of this variance). Hence it remains to answer how one should 
choose the vector a. The idea for this comes little bit clearer from the first 
order variation and Gorollarv 13.11 Indeed, as also pointed out in [22], the 
number D is the order of differentiability of X in the L 2 -sense. Hence if 
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D > 1, the variance must be at least of order (A tj) 2 so there is no hope 
to obtain even convergence in probability. Hence larger the D, larger the 
value of p should also be. Similarly, as 7 becomes close to number 2 it 
roughly means that D comes closer to 1 so once again one needs to refine 
the discretisation to obtain a Gaussian limit. More precisely, as 7 comes 
closer to 2 we see immediately that the variance is no longer enough to 

3 

compensate | 7 r n |2 in order to obtain central limit theorem. Hence one has 
to consider second order quadratic variations in order to stay in a Gaussian 
paradise. 


4 Examples 

This section is devoted to examples. We focus to reproduce some interesting 
and already studied examples to illustrate the power of our method rather 
than finding a complicated new examples, while now already it should be 
obvious how our approach can be used to study more complicated cases. 
More precisely, we study Brownian motion, fractional Brownian motion and 
related processes together with extensions sub-fractional Brownian motion 
and bifractional Brownian motion. Furthermore, our main focus is on first 
order quadratic variation for which we find sufficient conditions for the mesh 
to obtain almost sure convergence. In this context particularly interesting 
case for us is bifractional Brownian motion for which we are able to improve 
the sufficient condition proved in [28]. For simplicity we assume that the 
function k(n n ) is bounded. 

4.1 Standard Brownian motion 

Let X = W be a standard Brownian motion. Then it is known that the 
almost sure convergence holds provided that |7T n | = o (for recent re¬ 

sults on the topic see E3). Furthermore, this is a sharp in a sense that one 
can construct a sequence with | 7 r n | = O ^ 7 ^ 77 ^ such that almost sure con¬ 
vergence does not hold. Now the sufficiency of | 7 r n | = o ^ 7 ^ 77 ^ follows easily 
from concentration inequality (12.71) applied to the increments of Brownian 
motion. Indeed, in the case of standard Brownian motion the covariance 
matrix of the increments is diagonal, and we have 

||r ( ")||i = ||r (n) || 2 = |7T„|. 

Note also that if one uses Frobenius norm ||F^||f to obtain the upper 
bound, we have ||r^||^ = \/| 7 r n | provided that < C. Hence by using 

Frobenius norm we can only obtain half of the best possible rate even in the 
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case of standard Brownian motion. Finally, it is straightforward to obtain 
central limit theorem which, of course, is well-known already. 

4.2 Fractional Brownian motion 

Recall that a fractional Brownian motion B H with Hurst index H € (0,1) 
is a continuous centred Gaussian process with covariance function 

R(t, s) = i ( t 2H + s 2H -\t- s\ 2H ) . 

The case H = \ reduces to a standard Brownian motion. To obtain L p - 
convergence of general a-variations is straightforward by using Remark 12.21 

Proposition 4.1. Let B H be a fractional Brownian motion with H € (0, 1). 
Then there exists a constant Ch such that for a = jj we have 

E \ B S - B S-i i“ -* cht 


in L p for any p > 1. 

Remark 4.1. The exact value of the constant Ch is given by Ch = E|JV|h' , 
where N is a standard normal random variable. 


We now turn to the convergence of quadratic variation which is more 
interesting for us. Now it is natural to take <j>(x) = x 2 ^™ 1 , since for any 
partition of [0, T] we have 


E 

t^67T„ 


E (^-^_ x ) 2 
[A t k ]^-i 


= T. 


The following result is a direct consequence of Theorem 13.31 

Proposition 4.2. Let B H be a fractional Brownian motion with H G (0,1). 
Then 


(i) 


V B ■ = 

v n • 


7T„)-1 

E 

k =1 


(K - B, 


H ' 
tk — 1 ' 


[At fc ] 


2H—1 


in probability and in L p for any p > 1. Furthermore, the convergence 
holds almost surely for any sequence of partitions satisfying |7r„,| = 

° » Where 7 = maX (2 -=2Hi X ) ■ 


28 






(ii) There exists a constant C > 0 such that 


sup 


V B - T 

v n 


\JV ar(V B — T) 



P(Z < x ) 


< C|7T n | min (3’i- 2H ) 


where Z is a standard normal random variable. In particular, central 
limit theorem holds for all values H < |. 


Remark 4.2. Note that in the case H < \ we obtain similar sufficient 
condition as for standard Brownian motion. Indeed, the only difference is 
that since here the increments are not independent, we have to pose an 
additional assumption sup n>1 k(ir n ) < oo to obtain the optimal condition 

o (j- 5 — 

V log n 


Remark 4.3. By considering uniform partitions it can be shown that via 
concentration inequalities one cannot obtain any better result. It would be 
interesting to know whether the given conditions are optimal similarly as in 
the case of Brownian motion. However, for Brownian motion the counter¬ 
examples are constructed relying on independence of increments and, to the 
best of our knowledge, there exists no method to attack the problem for 
general Gaussian process. 


Remark 4.4. The limit theorems for quadratic variations of fractional 
Brownian motion are extensively studied in the literature. However, most 
of the related studies rely on uniform partitions and focus on generalisa¬ 
tions, e.g. to study Hermite variations or weighted variations rather than 
generalising the sequence of partitions. Furthermore, to recover the central 
limit theorem in the case H < | our approach is based only on simple linear 
algebra. For this reason our approach may be more applicable to generalise 
the results for arbitrary Gaussian processes while the obvious drawback is 
that it cannot provide a full answer to the problem. Finally, to the best 
of our knowledge the non-uniform partition are not widely studied in the 
literature. 


Remark 4.5. It is known that in the critical case H = f we also have 
convergence towards a normal random variable (see, e.g. [23] and references 
therein) with the only difference in the rate, i.e. the variance is of order 
instead of -k Given a priori knowledge that the variance is of order kfi-2: it 
is straightforward to recover also this critical case by Corollary 12.11 Hence 
again, it is sufficient to study the asymptotic behaviour of the variance. 
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4.3 sub-fractional Brownian motion 


The sub-fractional Brownian motion G H with parameter H € (0,1) is a 
centred Gaussian process with covariance function 

R(s, t ) = s 2H + t 2H -\ [s 2H + t 2H - \t - s\ 2H } . 

Note that as for fractional Brownian motion, value H = \ corresponds to a 
standard Brownian motion. 

Proposition 4.3. Let G H be a sub-fractional Brownian motion with H € 
(0,1). Then 


(i) 


U G := E 

67r n 


K - G L, 


T 


[A t k ] 2H ~i 

in probability and in L p for any p > 1. Furthermore, the convergence 
holds almost surely for any sequence satisfying |7r n | = o ^ (\ 0 g n )i ^, 

where 7 = max ^ 2 - 2 H ' ■*■) ■ 

(ii) There exists a constant C > 0 such that 

V?-T 


sup 

rrEM 


VMM 


< x — P(Z < X) 


< c |7r n r in (^§- 2i 0, 


where Z is a standard normal random variable. In particular, central 

3 

4’ 


limit theorem holds for all values H < 3 


Proof. Note that the case H = \ is already covered since it reduces back to 
the standard Brownian motion. Hence assume H 7 ^ -7 and let t > s. We 
have 


dstdcit, s) < C(t + s) 2H 2 + C\t-s\ 2H 2 


for some constant C and dc(s,t ) < C\t — s\ 2H . Now ( t + s) 2H ~ 2 < \t — 
s | 2 i /-2 f rom which the result follows immediately as in the case of fractional 
Brownian motion. □ 


Remark 4.6. We remark that the above result was already given in [2S] 
with the same rates although there the condition for the case H = ^ was 
| 7 r n || log | 7 r n || = which would follow from Remark 13.41 Obviously 

however, in this case we have standard Brownian motion so that | 7 r n | = 
o is sufficient. Note also that in this case one cannot obtain better via 

concentration inequalities. Indeed, this comes from the ’’fractional Brownian 
part” \t — s\ 2H . 
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4.4 Bifractional Brownian motion 


Particularly interesting case for us is Bifractional Brownian motion which 
was also studied in [28]. However, with our method we are able to improve 
the results of [28] , 

The bifractional Brownian motion is an extension of fractional Brownian 
motion first introduced by [21] and later analysed e.g. by [13] , Finally, 
bifractional Brownian motion was extended for values K G (1,2) such that 
HK G (0,1) by [I]. 

Definition 4.1. The bifractional Brownian motion is a centred Gaussian 
process B H,K with K = 0 and covariance function 

s) = ^ (( t 2H + s 2H ) K - I t- s \ 2HK ) 
with H G (0,1) and K G (0, 2) such that HK G (0,1). 

Remark 4.7. Note that K = 1 corresponds to ordinary fractional Brownian 
motion. It is straightforward to check that B H,K is H IF-self-similar and 
Hlder continuous of any order 7 < HK. For more details on the properties 
of bifractional Brownian motion we refer to [33] and references therein. 


While the main emphasis in [33] was integration via regularisation it was 
pointed out that one can prove that a-variations exists as a limit in L 1 in 
our sense. Hence the following result is obvious from remark 12.21 

Proposition 4.4. Let B H,K be a bifractional Brownian motion with H G 
(0,1) and K G (0,2) such that HK G (0,1). Then there exists a constant 
Ch,k such that for a = -jjj; we have 



- B 


H,K 
tk — 1 


->• Ch,kT 


in L p for any p > 1. 

Remark 4.8. In [33j the authors considered only the case K G (0,1]. How¬ 
ever, it is straightforward to obtain the claim for the case K > 1 by repeating 
the arguments. 

The next theorem studies the quadratic variation of bifractional Brown¬ 
ian motion. 

Proposition 4.5. Let B H,K be a bifractional Brownian motion with H G 
(0,1), K G (0, 2) and HK G (0,1). Then 
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(i) 


K 


H,K _ 


E 


B H ' K - B H,K 

tk—1 
2HK-1 


[A 4] 


2 1 ~ k T 


in probability and in L p for any p > 1. Furthermore, the convergence 
holds almost surely for any sequence of partitions satisfying |7r n | = 

(k^)> where ; 


• 7 = max ( 2 _ 2 HK , l) for K G (0,1]. 

* 1 = mm(l,2H)+l—2H K f 0r ^ ^ (■*■’ ^)- 

(ii) In the case K G (0,1] we have 


sup 

rrGM 



I—Kj-i 

- 2 1 ~ K T) 



— P(Z < x) 


< C\7T n \ min (^l- 2HK ) 


for some constant C and a standard normal random variable Z. In 
particular, central limit theorem holds for all values HI\ < |. 

(iii) In the case K G (1,2) we have 


sup 

xGM 


p 



1 —Kj~i 

- 2 l ~ K T) 



— P(Z < x ) 


< C\7r n \ min ^ 2H) ~ 2HK + ^ 


Proof. We assume K 1 since the case K = 1 reduces to ordinary fractional 
Brownian motion treated in Proposition 14.21 

The function d(s,t ) = P[B t — B s ] 2 is differentiable outside diagonal and 
we have 

d st d{s,t) = Ci|i - s\ 2HK ~ 2 + C 2 ( t2 H + s 2 H y- K 
for some unimportant constants C\ and C 2 . Furthermore, we have 

d(s, t) ~ \t — s\ 2HK 


as \t — s| —>• 0 which also corresponds to fractional Brownian motion. Now 
the term 1 1 — s\ 2HK ~ 2 can be treated as in Theorem 13.31 and for HK = tj 
this term vanishes. Consider next the term 

(t 2H + s 2H ) 2 ~ K ' 
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We are left to bound integrals 


h = 


f‘‘ 

J 0 Jtj—i 


(ts) 


2H-1 


x ( t 2 H + s 2Hy-K 


dsdt 


and 


ff 1 

v to v tj — 2 


/ 2 = 


(ts) 


2H-1 


rdsdt. 


,_ L (t 2 ^ -)- g2H^2—K 

We consider only /i since /2 can be treated similarly, and we denote by C 
any unimportant constant which may vary from line to line. By change of 
variables x = t 2H , y = s 2H and Tonelli’s theorem we have 


2 H H 


I\ = c 


/"T 

Jo Jt™ 


(xy) 


w( 2 ^- 1 ) 


—— X 2H 1 y2H 1 dxdy 


(x + y) 


2 H „f.2H 


= C 


r f 

Jo Jt‘ 


(x + y) K 2 d xdy 


j-i 
j. 2H n+ 2 H 


C [ [ (x + y) K 2 dxdy 

Jt? H , Jo 


= C 


Jo 


+y) K -'dxdy 


K- 1, 


3~ 1 


-C 


fT 


y* _1 dy 


Now for it > 1 we have y A 1 < C which leads to 


. ff-i 


dy < C(t 2H — < C|7r, 


imin(l,2H) 


7-1 


by the fact that for T > a > b > 0 and 7 € (0,1) we have a 7 — b 1 < (a — 6) 7 
and for 7 > 1 we have a 1 — W < C(a — b) by the mean value theorem. 
Similarly, for K < 1 we have 


y 


K -1 


1 ^ + 2HK + 2HK ^ 1 

dy < r 7 - - r 7 _i < |vr, 


imin(l,2i?ftT) 




jf -1 


Treating other integrals similarly the result follows by Theorem 13.31 with 
4>(x) = x 2HK ~ 1 . □ 

Remark 4.9. It may look like that for case K > 1 one gets better (i.e. 
larger exponent) by computing 




K -1 


dy = C(t 2 ^ A ' - g ™) < C|7T„r in(1 ’ 2 ^ X) . 
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t 2H 

However, this analysis cannot be used to cover, e.g. integral [ f 2 H (T + 
y) K ~ l dy. 

Remark 4.10. To compare our results with the existing literature, in [28] 
it was proved that almost sure convergence holds, in our notation, for value 
7 = for the whole range H £ (0,1) and K £ (0,1]. Note that by 

putting K = 1 and H = \ we have a standard Brownian motion, and this 
produces half of the best possible rate. Now in our result we have value 
7 = 2 - 2 HK w hfoh i s twice better compared to 1 _ 1 HK ■ Furthermore, we 
obtained even better for the range 2 HK > 1. Note also that, to the best 
of our knowledge, the case K > 1 is not studied in the literature before the 
present paper. 

Remark 4.11. Particularly interesting case is HK = In this case the 
quadratic variation exists in the ordinary sense which allows one to develop 
stochastic calculus with respect to B H,K although B H,K is not a semimartin¬ 
gale ES] if K £ (0,1). Now for this process we obtain similar condition to 
the one of standard Brownian motion. On the other hand, if K £ (1,2) 
and HK = 4, then the process B H,K is a semimartingale. However, in this 
case we only obtain condition \'K n \ 2H = o since K > 1 and HK = | 

implies H < i. 

Remark 4.12. In the case K £ (0,1] we obtain a sufficient condition HK < 
~ for the central limit theorem to hold which of course is not at all surprising. 
Similarly, in the case K £ (1,2) we obtain sufficient condition HK < | 
provided that H > \. However, in the case H < \ something odd seems 
to happen. Indeed, if HK >1 then 2 H + -( — 2 HK < 0 so that the given 
Berry-Esseen bound does not converge to zero. On the other hand, even if 
HK < | it is not necessarily true that 2 H + tj — 2 HK > 0 so that condition 
HK < | is no longer sufficient. Indeed, even in the semimartingale case 
2HK = 1 we have 2 H + \ — 2HK < 0 for values H £ (0, |]. 

Acknowledgements Lauri Viitasaari was partially funded by Emil Aal¬ 
tonen Foundation. 


A Proof of Lemma 12.2 


A simple application of (12.41) yields 


EV n = E ( E 

k,j =1 


k j 
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Next we compute EV^. We have 

n 

< = e n 



(A.l) 


Recall next that all information of a Gaussian vector is encoded to the 
covariance matrix r^ n ) so that k- moments of a centred Gaussian vector 
(Y \, Y' 2 ,..., Y n ) can be computed via formula 

E[Y^Y^ . ..Y^] = ^E [Y^Y^] .. .E K^Y^}, 

a 

where the summation is over all permutations a of numbers {1,2,... ,n}, 
hence producing n! terms. Applying this to vector 8-dinrensional vector 
(yl n) , y / J n) , Y- n ^...., Y l (n) 'j and taking expectation on 



2 





2 


we obtain terms of form 


Ai(cj) = 

A 2 (a) = E 
A 3 (cr) = E 
A^(a) = E 

A 5 (cr) = E 


E 


Y 


\r( n ) \r{ n ) 
Y cr(l) Y a{2) 


2 r 


E 


\r( n ) \ r ( n ) 

M3)-V(4) 


1 2 


(n) 

(1) 


1 2 


E 


Y 


(n) 

( 2 ) 


1 2 


E 


\s( n ) \ r ( n ) 

M3)G(4) 


11 2 


Mi) 

2 

E 

M2) 

2 

E 

\s( n ) 

tr(3) 

2 

E 

\M) 

M 4) 


1 2 


w( n ) yd n ) 

E 

xA n ) y( n ) 

E 

yA n ) yd n ) 

E 

y( n ) 

. cr(!) cr ( 2 )_ 


cr(2) cr(3) 


cr(3) ct(1) 


<r( 4 )_ 


1 2 


y( n ) y( n ) 

E 

y( n ) y( n ) 

E 

-1 

E 

y( n ) y( n ) 

cr(l) <r(2) 


cr(2) cr(3) 


L o-(3j o-(4)j 


cr(4)' r fr(l) 


and 


where a = (cr(l), <r(2), cr(3), cr(4)) can be any permutation of indices 
{i,j, k, l}. Next we note by symmetry of covariance and summing over 
symmetric set {1 < i,j, k,l < n} we obtain that for each p = 1,..., 5 and 
any permutation a we have \( a ) = Yli,j,k,i= l A p (cto), where cr 0 is 

any fixed permutation. For example, we obviously have 


n 

E [ E 

y( n )y( n ) 

* J 

- 2 

E 

y( n )y( n ) 

O n 

] = E [ E 

IT 

i_i 

- 2 

E 

i 1 


i,j,k,l =1 
Consequently we obtain 
2 


i,j,k,l=l 


E E 

i,j,k y l=l 


Y 


(n) 


Y 


(n) 


(n) 


Y, 


(n) 


y ] ^ ] a pA P (c jq) 
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for arbitrary reference permutation ao and some weights a = (ai,...,a 5 ). 
Note also that the weights a p ,p = 1,..., 5 are independent of indices i. j, k, l 
and the underlying Gaussian process. Now treating rest of the terms in 


n 


pe{i,j,k,l} 


Yr 


(n) 


-E Y r 


r(n) 


similarly we conclude that 


ek 4 = 


5 


}^ b P A p( a o) 

i,j,k,l=1 p= 1 


with some weights b = (b\, ..., 65 ) independent of i,j, k , l and the underly¬ 
ing Gaussian process. Next we claim that b = (12,0,0,0,24). Of course 
the weight vector b could be computed exactly via combinatorial arguments 
but we wish to use a more subtle argument by relying on the classical cen¬ 
tral limit theorem for a sequence of independent standard normal random 
variables. We begin by computing the values 64 and 65 which are rela¬ 
tively easy to compute directly. First note that terms A$ are produced only 


by the term E 


Y, 


(n) 


Y 


(n) 


(n) 


Y, 


(n) 


Furthermore, terms of 


form A$ are produced by permutations of indices {i,j,k,l} which gives 
65 = 4! = 24. Consider next 64 . Terms of form A 4 are produced from 


E 


Y, 


(n) 


Y 


(n) 


Y 


(n) 


(W)' 


by first picking one variable, Y, 


(n) 


say, to get E 


Y, 


(n) 


and then organising the remaining three into 3! = 6 


ways which produces 4! = 24 (the first one can be picked in 4 ways). On 
the other hand, computing the product (1A.1D we obtain 4 terms of form 


E 


Y, 


(n) 


Y 


(n) 


(n) 


Y 


(n) 


and with similar analysis we obtain 


that each term produces A 4 exactly 3! = 6 times. Due to the minus sign 


in terms —E 


Y 


(n) 


and the fact 24 — 4 x 6 = 0 we obtain 64 = 0. It 

(n) 

remains to prove that b± = 12 and 62 = 63 = 0. For this purpose let Y p ' be 
a sequence of independent standard normal random variables Yfc. Then by 
the classical central limit theorem we have 


S n : = 


1 


n 






in distribution and consequently, ES 1 ^ ->-3. In this case we have 

n 

bi ^2 Ai(a 0 ) = hn 2 

i 


n 


h 2^ ^3(co) = b 2 n 3 
i,j,k,l=l 
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and 


h ^2 ^Oo) = b 3 n 4 

i,j,k,l =1 

so that 

ES^ = —^ [bi™ 2 + fan 3 + b 3 n 4 ] . 

Now since bk,k = 1,2,3 is independent of n and the underlying Gaussian 
process, the convergence ES^ —>• 3 implies b 2 = b 3 = 0 and b\ = 12. This 
completes the proof. 
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